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special relativity are found for four geometrically defined relative velocities in 
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1 Introduction 


In special relativity, the velocity addition formula for three inertial observers 
whose relative motions are spatially collinear may be expressed as, 


Vi + V2 
1 + V1V2 ' 


( 1 ) 


Here and below, relative to a central observer, vi is the velocity of a secondary 
observer; V 2 is the velocity of a test particle relative to this secondary observer; 
and V 3 is the velocity of the test particle relative to the central observer. An 
alternative and sometimes useful equivalent relationship is. 


1 - 'Ca _ / 1 - ni \ / 1 - U2 \ 
1 + ^3 Vl + ui/ Vl + ^ 2 / 


( 2 ) 


Are meaningful velocity addition formulas, analogous to these, possible in gen¬ 
eral relativity? At first glance, the answer appears to be no. General relativity 
provides no a priori definition of relative velocity when the test particles and 
observers are located at different space-time points. Different coordinate charts 
give rise to different conceptions of relative velocity. 


To avoid such ambiguities, we study the coordinate independent, purely ge¬ 
ometrically defined Fermi, kinematic, spectroscopic, and astrometric, relative 
velocities introduced by V. J. Bolds in [T], and subsequently developed in [5], 
0, la, 0, i, la, 0, M- The definitions of these geometric relative velocities 
depend on two distinct notions of simultaneity: “spacelike simultaneity” (or 
“Fermi simultaneity”) and “lightlike simultaneity.” 


The Fermi and kinematic relative velocities depend on spacelike simultaneity. 
Two events are simultaneous in this sense if they lie on the same Fermi space 
slice Mt determined by a fixed Fermi time coordinate r. To define Fermi 
coordinates (see, e.g., m), consider a foliation of some neighborhood U (which 
might be the entire spacetime) of a central observer’s geodesic worldline, Po{t), 
by disjoint Fermi spaceslices {Mr} defined by, 

Mr = ^;H0)- (3) 

Here the function y;,- : K is given by, 

V’rip) = g{exp~^^^p, $o{t)), ( 4 ) 

where the overdot represents differentiation with respect to proper time t along 
/3o, g is the metric tensor, and the exponential map, expp(z;) denotes the evalu¬ 
ation at affine parameter 1 of the geodesic starting at point p G M, with initial 
derivative v. In other words, the Fermi spaceslice Mr of all r-simultaneous 
points consists of all the spacelike geodesics orthogonal to the path of the cen¬ 
tral observer /3o at fixed proper time t . A convenient coordinate system on Mr 
consists of two angular coordinates together with the proper distance p from 
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the central observer’s path to a point on A4r (see Hi). 

The Fermi relative velocity, Fpermi) of a test particle relative to an observer, 
(io{T), is a vector field along orthogonal to the 4-velocity U of the ob¬ 
server (and therefore tangent to at each proper time r of the observer. 

For a test particle undergoing purely radial motion (which is our focus), the 
magnitude of Fpermi, or the Fermi relative speed, upermi, is the rate of change of 
proper distance p of the test particle away from the observer / 3 ( t ), with respect 
to proper time r. Eq. 0 is thus a statement about Fermi relative velocities in 
the context of special relativity. For more general motion (i.e. non radial), the 
definition of Fermi relative velocity may be found in [T] . 

The kinematic velocity of a test particle at a spacetime point g, relative to the 
central observer /3o(t), is found by first parallel transporting the test particle’s 
4-velocity u’ along a radial spacelike geodesic (lying on a Fermi space slice) to 
the 4-velocity denoted by Tqpu' in the tangent space of the central observer at 
spacetime point p = /3o(to). If the 4-velocity of /?o at the point p is u, then the 
kinematic relative velocity Ukin of the test particle is defined to be the unique 
vector orthogonal to u, in the tangent space at p, satisfying Tqpu' = ^{u + Ukin), 
where the scalar 7 is uniquely determined and is given by 7 = 1 / 

The spectroscopic (or barycentric) and astrometric relative velocities can be 
found, in principle, from spectroscopic and astronomical observations. These 
two relative velocities depend on lightlike simultaneity, according to which two 
events are simultaneous if they both lie in the past-pointing horismos at the 
spacetime point p of the central observer (the past-pointing horismos is tangent 
to the backward light cone). More specifically, using the notation above, let 
i/'t : W K be given by, 

V't (p) = g(exp“ p, p), (5) 

where in this context, we denote proper time along /3o by T0 Now define 

£;, = ^-1(0)-{/3o(t)}. (6) 

The 3-dimensional submanifold, Et, is called the horismos submanifold at the 
spacetime point /3o (t) (where proper time T is fixed). An event q is in Er if and 
only if g /3o(t) and there exists a lightlike geodesic joining /3o(t) and g. 
has two connected components, Ejf and E~, respectively, the future-pointing 
and past-pointing horismos submanifolds of (3o{t) (see [TT] and the references 

^In Sect|^ the symbol T will denote the time coordinate in optical coordinates which is 
different from the Fermi time coordinate r. However, these two coordinate times are identical 
along the path /do where they are both proper time. 
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therein). 

We note that for a comoving observer /3o, it was proved in [6] (see also [12]) 
that for a broad class of Robertson-Walker cosmologies, the maximal chart for 
Fermi coordinates is exactly the causal past of /3o so that, 

IJm, = 1Jf;-u/3o. (7) 

r>0 T>0 

The spectroscopic relative velocity r^spec is calculated analogously to Ukin, de¬ 
scribed above, except that the 4-velocity u' of the test particle is parallel trans¬ 
ported to the tangent space of the observer along a null geodesic lying on the 
past-pointing horismos of the central observer, instead of along the Fermi space 
slice. The astrometric relative velocity, Uast, of a test particle whose motion is 
purely radial is calculated analogously to upermi) as the rate of change of the 
affine distance (see Sect. [^, which corresponds to the observed proper distance 
(through light signals at the time of observation) with respect to the proper 
time of the observer, as may be done via parallax measurements. A complete 
description is given in [T]. 


In Minkowski spacetime, the coordinates of three Lorentz frames (with relative 
velocities as in Eq. 0 ) in “standard configuration” relative to each other are 
related by Lorentz boosts along a given space axis of a central observer, and 
all origins of coordinates coincide when all of the time coordinates equal zero. 
The natural generalization to Robertson-Walker cosmologies is to consider three 
observers co-moving with the Hubble flow, with collinear relative motion (i.e. 
with two fixed space coordinates). In the case of the Milne universe, this config¬ 
uration amounts to three special relativistic observers in standard configuration. 


In this paper, we obtain generalizations of Eq Q for comoving observers in 
Robertson-Walker spacetimes. As a special case, using general methods we re¬ 
cover Eq Q for the Milne universe which is diffeomorphic to the forward light 
cone of Minkowski spacetime under a simple change of coordinates. In the Milne 
universe, Eq. 0 holds for all of the geometrically defined relative velocities ex¬ 
cept for the astrometric velocity, which follows a different addition law (Eq. 
(52)). The addition laws for the four geometrically defined relative velocities 
are strikingly different. For example, we show that for lightlike simultaneous 
observers, the spectroscopic relative velocity follows Eq ([^ in all Robertson- 
Walker spacetimes considered in this paper, while the addition formulas for the 
other relative velocities are quite different. 


Unlike Minkowski spacetime, for a general Robertson-Walker cosmology, the 
distinct foliations by sets of simultaneous events for each of the three observers 
forces the existence of more than one velocity addition formula, depending on 
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how comparisons are made. 


This paper is organized as follows. Sections anddevelop notation for Fermi 
and optical coordinates and spacetime paths of observers and test particles. In 
Sections and we express the Robertson-Walker metric in Fermi coordinates 
and optical coordinates and slightly generalize formulas for relative velocities 
found in |3] that we use in the proofs of our theorems. Sections and give 
general velocity addition formulas for the Fermi, kinematic, spectroscopic, and 
astrometric relative velocities. In Section]^ we show how, in special cases, one 
may use relationships between the Hubble velocity and the four geometrically 
defined relative velocities to derive velocity addition formulas. We use that 
to derive velocity addition formulas for Robertson-Walker cosmologies with any 
power law scale factor. Examples for the de Sitter universe and particular power 
law cosmologies are given in Sectionand Section]^ gives concluding remarks. 


2 The Robertson-Walker metric 

The Robertson-Walker metric on space-time A4 is given by the line element, 

ds^ = -df + a?{t) [dx^ + Sl{x)dVt^] , (8) 

where -|- sin^ 9 d(p^, a(t) is the scale factor, and, 

{ sin X if fc = 1 

X if fc = 0 (9) 

sinhy if fc = — 1. 

The coordinate t > 0 is cosmological time and X) d, (p are dimensionless. The 
values -|-1,0, —1 of the parameter k distinguish the three possible maximally 
symmetric space slices for constant values of t with positive, zero, and negative 
curvatures respectively. The radial coordinate x takes all positive values for 
fc = 0 or —1, but is bounded above by tt for fc = -1-1. 

We assume throughout that A: = 0 or —1 so that the range of x is unrestricted. 
The techniques employed for these two cases may be extended to the case A: = -|-1 
with the additional restriction that x < tJ" so that spacelike geodesics do not in¬ 
tersect. 

There is a coordinate singularity in ([^ at x = 0, but this will not affect the 
calculations that follow. Since our intention is to study radial motion with 
respect to a central observer, it suffices to consider the 2-dimensional Robertson- 
Walker metric given by 

ds2 =+ a2(^)dx^ (10) 
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for which there is no singularity at y = 0 and we may allow y £ (— 00 , 00 ). We 
assume throughout that a{t) is a smooth, increasing function of f > 0. 


3 Notation for three observers 

We denote a central observer by the path f3o{to) = ito,0), a secondary observer 
by /3i(ti) = (ti,Xi)) and a test particle by /32{t2) = {t 2 ,X 2 )- Each spacetime 
path is parameterized by its proper time, and we take /3i and /32 to be comoving 
so that xi and X 2 are constant. Without loss of generality we take X 2 to be 
positive, while allowing xi to assume both positive and negative values. We 
allow for the possibility that xi > X 2 - 


The Fermi, kinematic, spectroscopic and astrometric velocities of a test particle 
relative to a particular observer, /3(t), are smooth vector fields defined along 
the path /3(t). Those vector fields (along paths of observers) are denoted by 
upper case letters, and scalar fields (relative speeds) by lower case letters. 


In what follows, the relative velocity vector fields are defined on the spacetime 
paths j3{j and /3i, i.e. the central and secondary observers. From their definitions 
(see [T]) all four relative velocity vector fields are spacelike and orthogonal to 
the 4-velocities of the observers. They are each multiples of the unit vector field 


a{t) dx 


( 11 ) 


The Fermi, kinematic, spectroscopic and astrometric velocity vector fields are 
denoted respectively by Vpermi = ?^Fermi5, I4in = 'ykin5,I4pec = VspecS, and 
East = VsstS. We note that the scalars t^Fermi, ^^kin, etc. can be positive or 
negative, indicating velocities in the same or opposition direction as S. We note 
that this notation differs from that used in [J]. 


Subscripts will be used to distinguish the different relative velocities. Subscripts 
1 , 2 and 3 will denote a velocity of j3x relative to /3o, h relative to /3i and P 2 
relative to /3o, respectively. E.g. WFermi 2 is the Fermi velocity of the test particle 
relative to the secondary observer. 


4 Fermi coordinates and spacelike simultaneity 


It may be shown that the metric (101 expressed in Fermi coordinates (r, p) for 
the central observer, /3o, (whose y-coordinate is taken to be zero) has the form, 


ds^ = 5rr (f, p)dr^ + dp^. 


( 12 ) 
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General formulas for g^-r were derived in uni. Non negative values of the 
spatial coordinate p give the proper distance along spacelike geodesics orthog¬ 
onal to the world line of the comoving Fermi or central observer. However, 
consistent with Eq. (101 we also allow p to take negative values. 


Formulas for relative velocities developed in nil] can be easily extended to 
accommodate negative values of the spatial coordinates. At a given proper time 
r of the central observer, suppose that the Fermi coordinates of a test particle 
are (r, p) (and the coordinates of the central observer are /3 o(t) = (''"jO))- Let 
the curvature coordinates of the test particle be (t, x) (where y is fixed). Then 
the kinematic velocity of a comoving test particle relative to the central observer 
is given by. 


Ukin = sgn(x) 



a2(r)' 


(13) 


The Fermi speed is given by upermi = dp/dr for a radially moving test particle. 


The common dependence on spacelike simultaneity of the Fermi and kinematic 
relative velocities allows for a direct comparison of these two notions of relative 
velocity of a test particle at a given spacetime point. The Fermi and kinematic 
speeds of a radially moving test particle at the spacetime point (r, p) in the 
Fermi coordinates of a central observer are related by, 

'i^Fermi — \/ Prr (f, p) Upin. (1^) 

In what follows we will need Fermi coordinates for a comoving observer with 
some fixed y-coordinate, y^-, not necessarily equal to zero. As in the previous 
section, let j3j represent the worldline of a comoving observer whose y coordinate 
is Xj ■ The proper time of /3j is again Tj , and Fermi coordinates for the observer 
Pj may be constructed. In this case we denote the leading metric coefficient by 
gTjTj, so that, for example, may also be written as Proxo- 

Functional relationships among the coordinates, Tj, yj, t, and y, follow from 
the observation that the vector field. 


X = -sgn(y - Xj) 



_ a a{Tj) d 
dt a?{t) 9y’ 


(15) 


is geodesic, unit, spacelike and orthogonal to the 4-velocity of /3j, i.e. X is 
tangent to for each value of Tj. Let q = {t, y) G Xtrj ■ Then there exists an 
integral curve (for the vector field X or —X) from p G Pj to p, and thus using 
Eq. (15) we can find a relationship between Tj, yj, t, and y (see [4]): 
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1 


(16) 


a{Tj) 


i{t) 


a^(Tj) — a^(t) 


-.di= \x-Xj\- 


Remark 1. Eq.{13\ can be easily adapted to yield velocities of a test particle 


relative to the comoving observer, fdj by replacing sgn(x) by sgn(x — Xj) o-nd r 
byTj. 


5 Optical coordinates and lightlike simultaneity 

In the framework of lightlike simultaneity, it will be convenient to use optical 
coordinates (also known as observational coordinates). Following [3], we may 
express the Robertson-Walker metric in the optical coordinates, (t,(5), deter¬ 
mined by the central observer /3o, 


ds^ = 5TTdT'^-|-2sgn((5) dxd^ = —2 1-rwl<5| ~ x 


i(T)| 


1 of (t) 


z(t) 2a2(T) 


where t(T, 6) is given implicitly by, 


dT^-|-2sgn(i5) dT dd, 

(17) 


<5 = <5(i,x) = sgn(x)^ 


T(t.x) 


l(u) 


It a{T{t,x)) 
and T{t,x) is determined implicitly from the equation, 


du. 


1 


i{u) 


du = Ixl- 


(18) 


(19) 


Remark 2. Along the path, fdo, of a comoving central observer, the three time 
coordinates t,r, T (respectively from curvature coordinates, Fermi eoordinates, 
and optical coordinates) take the same values and are identieal. They differ, 
however, at spacetime points off the path /Sq. We note that our notation for the 
optical time coordinate T differs from that used in We have also dropped 
the subscript tint and x used in that reference. 

Formulas for relative velocities developed in [3] can be readily extended to ac¬ 
commodate negative values of the spatial coordinates. The spectroscopic and 
astrometric velocities of a comoving test particle at fixed spatial coordinate x 
relative to the central observer, /3o(t), are given by, 


Vspec = sgn(x) 


a^(T) — a^{t) 
a^(T) + a?{f) 


and 


Uast = sgn(x) 1 - |5| 


i(t) af{t) 


i(t) a2(T) 


( 20 ) 


( 21 ) 












where 6 is the affine distance parameter from the central observer to the co¬ 
moving test particle. As in the previous section, the formulas here for relative 
velocities can be easily adapted to yield velocities of a test particle relative to 
the comoving observer, Pj by replacing sgn(x) by sgn(x — Xj) and r by Tj. 

The following comparison of the astrometric and spectroscopic relative veloci¬ 
ties is also possible (see i) because of their common dependence on timelike 
simultaneity. 


_ sgn((5) 

"^ast — 


6) + 


1 - Sgn(^) Uspec \ 
1 -f Sgn(5) VspecJ 


( 22 ) 


6 Velocity addition formulas for spacelike simul¬ 
taneity 

In this section we derive velocity addition formulas for the Fermi and kinematic 
relative velocities of comoving observers. In general there can be no single for¬ 
mula like Eq. Q for either relative velocity. This is because the respective Fermi 
submanifolds of simultaneous events for the central and secondary observers do 
not coincide, and as a consequence the spacetime points at which relative veloc¬ 
ities are calculated can be reasonably chosen in more than one way. We show, 
however, for the special case of the Milne universe (essentially for the case of 
special relativity), that Eq.Q follows as a special case of Theorem]^ and Corol¬ 
lary below. 

The figures below depict three general scenarios for which addition formulas 
may be deduced, and we refer to them in Theorem and Corollary In 
each hgure, the comoving paths Pi as described in Sect, [^are indicated. The 
vertical axis is cosmological time, which is also proper time for each comoving 
path, and the horizontal axis is the x-axis. The curves labelled by 'kg or dtj 
are spacelike geodesics orthogonal to the central observer Pq; tki is a spacelike 
geodesic orthogonal to the secondary observer Pi. 
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Figure 1: Scenario I. Here 'I'o is the spacelike geodesic orthogonal to /3o, and 
is the spacelike geodesic orthogonal to /3i. is the unique event tko n /3i 

and (t^,X 2 ) is the unique event tpi n /32- The velocities of j32 and /3i relative 
to /3o are taken at proper time Tq of /3o, while the velocity of /?2 relative to (3i is 
taken at proper time ti of /3i. Note that in this case the velocity of /?2 relative 
to the observers will be computed for separate events on the world line of ^ 2 , 
i.e. the time discrepancy has been placed on the world line of the test particle. 



Figure 2: Scenario II. Here and are the spacelike geodesics orthogonal 
to /3o at f = Tq and t = Tq , respectively. ( ti , xi) is the unique event dtj^ n I3[, 
and d^i is the spacelike geodesic orthogonal to /3i emanating from this event. 
Pi, d^Q and dti intersect at the event (t 2 ,X 2 )- Note that in this case the time 
discrepancy is on the world line of the central observer. 
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Figure 3: Scenario III. Here 4'o is a spacelike geodesic orthogonal to /3o, and 
intersects both j3i and /32 at (r^^, xi) and (t 2 , X 2 ), respectively, is the unique 
geodesic orthogonal to /3i such that (t 2 , X 2 ) = /32 H 'I'l. In this scenario the time 
discrepancy is on the world line of the secondary observer. 


Remark 3. Figures\^ a 7 T,(i[^ depict the case that X 2 > Xi > 0. However it 
is possible to connect the spacetime events shown in each diagram by spacelike 
geodesics for arbitrary configurations of the observers and test particle. The 
results we derive are thus valid in greater generality than the figures imply. 

Theorem 1. For the scenarios depicted in Fi^rtres and let Ukini be the 
kinematic speed of pi relative to Po; nkin 2 be the kinematic speed of p 2 relative 
to Pi; and Ukins be the kinematic speed of P 2 relative to Po. Then, 


(1 ^kins) ~ ^ 2 !^+! ^kinl)(l ^kin2)i (^3) 

where the proper times t'^ and t~ are indicated in the figures and are determined 
by the indicated spacelike geodesics. 


Proof. For fig. using Eq. (13), we may write, 


^kinl — 1 ~ 


'M. 

Kro) 


^kin2 — 1 ~ 


^(n) 


■^kinS — 1 ~ 


\t-) 

^{ro) 


(24) 


The result now follows by direct substitution of Eqs.(24| into Eq. 
for figures and we may write 


Similarly 


^2 _ ^ _ aHn) 

“ a^t-) 


^kin2 


= 1 - 


a2(Ti) 


^kin3 


= 1 - 


a^(T2) 

a 2 (t+) 


(25) 
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and 


'^kinl 


= 1 - 




a^(T2) 


'^ki„2 - 1 ^2(i- 


^kin3 


= 1 - 


a^ij2) 

a^(T-o)’ 


(26) 


respectively. Substitution of Eq. (251 or Eq. (261 into Eq. (23) verifies the 
identity. 

□ 


Using Eq.(14|, the following corollary is immediate. 


Corollary 1. Following the notational conventions of Sect^ the Fermi relative 
velocity addition formula for comoving observers is given by, 



'^FermiS 

9toto (s) 


) f. , ^lermil \ f. , ^iermi2 \ 

I 9roro{q)) \ 9r,rAs)J’ 


(27) 


where s and q stand for the appropriate Fermi coordinates for the spacetime 
points indicated in Figures^ and where grin (s) is replaced by grin {s') 

for the scenario depicted in Figure Q 

Example 1. The velocity addition formula for special relativity (more precisely 
for the Milne Universe) is a special case of Theorem^ We illustrate this for 
the configuration shown in Similar calculations show that the result holds 

for the configurations shown in Figs. [7] and The transformation formulas 
from curvature coordinates {t, x) to Fermi coordinates (r, p) are, 


T = t cosh X P = t sinh y 


(28) 


(see e.g. m)- Expressed in Fermi coordinates, the metric for the Milne Universe 
with k = —1 and aft) = t) becomes the Minkowski metric, i.e., in two 
spacetime dimensions, 

ds^ = —dr^ 

Applying Eq. (28) to gives, 


and 


Since a{t) = t, we then get. 


i{t+) 


cosh(y2) 


+ dp^. 

(29) 

= T 2 cosh y 2 

(30) 

12 - Xl)- 

(31) 

= 1 + tanhyi tanh(x 2 - Xi)- 

(32) 


a(ti ) cosh(y 2 - Xi) coshyi 
By combining Eqs. (30) and © with Eq.^^, this may be expressed as. 
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(33) 


aHtt) 

a2(t-) 


= (1 + Wkinll’kin2) 


2 


Now substituting Eq. 


into Eq. (23) gives, 


(1 + Wki„iWki„2)^ (1 - vLa) = (1 - '(^Li)(l - «L2), (34) 


which reduces to 


Vkin3 = 


^kinl + 'i^kin2 
1 + l’kinll’kin2 


(35) 


Applying Corollary^ with grr = — 1 for the Milne universe then recovers 


'I’Fermil + 'f^Fermi2 

^FermiS — 

i + 'l’Fermil‘^’Fermi2 


(36) 


We note that for a comoving observer with fixed coordinate x, it follows directly 

VFermi = tanhx- (37) 


from Eq.{28) that, 


Thus, X is the rapidity parameter for the Lorentz boost determined by VFermi (see, 
e.9; m), and this Lorentz transformation is a hyperbolic rotation by angle X- 
In it was also shown that Wkin = "yFermi = tanh x for a comoving test particle 
in the Milne universe, or, in the language of special relativity, a test particle 
with constant velocity departing from the central observer at time zero. 


7 Velocity addition formulas for lightlike simul¬ 
taneity 

In this section we derive velocity addition formulas for the spectroscopic and as¬ 
trometric relative velocities of comoving observers. For the observers /3o,/3i,/32, 
the prototype configuration is depicted in Fig|^ in which the spacetime events 
p, q, and s are all lightlike simultaneous. 
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Figure 4: Elements involved in the study of velocity addition for simultaneous 
(in optical coordinates) comoving observers. A is a lightlike geodesic. 


Theorem 2. For simultaneous (in optical coordinates) comoving observers the 
spectroscopic relative velocities are related by the special relativity addition for¬ 
mula, Q, i.e., 

^specl 4“ 'Cspec2 

^^spec3 = T—- (38) 

^ I ^specl^spec2 

Proof. For simultaneous comoving observers the spectroscopic velocities de- 
picted in Figare given by (20) as, 


^specl — 


a'^ito) - a^(ti) 




a 2 (ti)’ 


^spec2 — 


a^(ti) - a^(t2) 
a^ih) + a^h)’ 


and 


^spec3 — 


a^(to) - a^{t2) 
a'^{to) + a?{t2) 


(39) 


(40) 


The result follows from direct substitution of Eqs. (39) into the right side of 
Eq. (p^ which then becomes Eq. (401. □ 


Theorem and Fig|^ assume that X 2 > Xi > 0- However, this restriction is not 
essential. For any ordering of the spatial coordinates of the observers and test 
particle a measurement scheme can be found so that the spectroscopic velocities 
are related by Eq. ([^, but the measurements may not all be simultaneous. To 
illustrate, we prove this for the case that xi < 0 < X 2 - Fig. [^indicates how the 
spectroscopic relative velocities are to be measured for this case. 
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Figure 5: In this scenario, UspecS and nspec 2 are the spectroscopic velocities of 
j 32 at the spacetime point s relative to /3o and /3i, respectively, rasped is the 
spectroscopic velocity of /3i at the spacetime point q relative to /3o. Ai and A 2 
are lightlike geodesics. 


Theorem 3. The spectroscopic relative velocities of the comoving observers 
depicted in Figure are related by the special relativity addition formula, Eq. 
Q, i.e., 

^^pec3 ^ 

-L H" '^specl'^spec2 


Proof. For this scenario t'spec 2 and UspecS are given in Eqs. 
spectively, while Uspeci is given by 


(39) and (40), re- 


/ a^(ti) - a^(to) \ _ a^(to) - a^(ti) 

\a'^{to) + a'^{ti) J a2(to) + a2(fi) ■ 


Thus the spectroscopic velocities are given by the same expressions as those 
used in the proof of Theorem and the result follows exactly as before. □ 


Using the velocity addition formula for the spectroscopic velocity and the func¬ 
tional relationship between the astrometric and spectroscopic velocities given 
by Eq. (21), we derive an addition formula for astrometric velocities for the 
lightlike simultaneity scenario depicted in Fig. 


Lemma 1. Suppose that the spacetime points Polfo), /3i(ti), and ^ 2 ( 12 ) He on a 
null geodesic in the past pointing horismos Eff. Denote the optical coordinates 
of j3i{ti), and / 32 (t 2 ) relative to the central observer (do emd (to:^ 2 ) 

respectively, and assume 82 > > Q as in Fig. Let 621 be the affine distance 

coordinate of P 2 {t 2 ) in the optical coordinate system for fi (i.e. regarding fi as 
the central observer). Then, 


62 


Si + 


aifi) 

a{to) 


821 


(43) 
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Proof. From Eq. (18 1 with x > Oj 


<52 = 


1 


a{to) 

1 

a{to) 


rto 


't2 


= Si 


ajti) 

a(to) 


r(u)dr 


r(u)dr 


J21- 


a(^j_ 

a(to) a{ti) 


i(M)dr 


(44) 

(45) 

(46) 

□ 


Theorem 4. Following the notation o/Femma [7| the astrometric relative ve¬ 
locities of the simultaneous comoving observers depicted in Figure^ are related 
by 

— (2Wast3 + ffTT(^07 <^ 2 )) = (Sriastl + 5tt(^0j l^l)) (2fast2 + 5tiTi (^li <52l)) ; (47) 


where ^tiTi is the leading metric coefficient in optical coordinates for the sec¬ 
ondary observer fi and 

<521 = ^(<52-^1). (48) 

a{ti) 

Proof. Taking into account Remark]^ and combining Eqs. (22) and (38) (see 
also Eq. ©) gives, 


■Casta — ( ffTT(io7^2) + 


1 - U 


sped 


1 - u 


spec2 


It follows from Eq. (221 that, 

1 '^sped 


and 


1 't^sped 
1 '^^spec2 


1 H“ ‘^^sped 1 H“ ‘?^spec2 


I-\-u, 


spec2 


— —2Uast2 — 5 tiTi (il7 <52 i) • 


The result now follows by combining Eqs (49), (50), and (51). 


(49) 

(50) 

(51) 

□ 


Example 2. For the Milne Universe, it is easily verified from Eq. ( |17| ) that 
^tt = ~1- Applying Eq. (47) then gives the following special relativistic veloc¬ 
ity addition formula for the astrometric relative velocity, for the configuration 
depicted in 

ast^ — 'Casti 4“ ast2 ast\'^ast2 • (^2) 

IFe note, however, that different velocity addition formulas for the astrometric 
relative velocity hold for configurations different from that depicted in F^g\^ 
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8 Hubble Additivity 

In this section, we consider Robertson-Walker spacetimes for which the geomet¬ 
rically defined velocities of comoving test particles can be expressed as functions 
of the Hubble velocity. Examples of such spacetimes are given in the following 
section. In this situation, the addition formulas for the geometrically defined 
velocities may take simpler forms than for more general scenarios. 

The usual foliation of a Robertson-Walker spacetime by maximally symmetric 
space slices {Et}, parameterized by cosmological time t, determines a notion 
of simultaneity different from spacelike or lightlike simultaneity and naturally 
leads to the Hubble velocity and Hubble’s law, 

V = d{t) = a{t)x = Hd. (53) 

Here v is Hubble velocity, H is the Hubble parameter, d = a{t)x is the proper 
distance on Ej from the observer to the comoving test particle with coordinate 
X, and the overdot signifies differentiation with respect to t. 


The velocity addition formula for the Hubble velocity of comoving objects, rela¬ 
tive to a comoving observer is particularly simple. Following the notation of the 
previous sections, let the secondary oberver, /3i, whose hxed spatial coordinate 
is xij have Hubble speed vi relative to /3o- The Hubble speed of /32, at X 2 , 
relative to /3o is denoted by V 3 , and V 2 denotes the Hubble speed of relative 
to /3i. Then at cosmological time t. 


W 3 = d(r)x2 = d(T)[xi + (X2 - Xi)] = vi + V 2 . (54) 

We see that the velocity addition formula for the Hubble velocity is Galilean. 


In general (see [3]), all four of the geometrically defined relative velocities of a 
comoving test particle are uniquely determined by the two coordinates, (t, x), 
where x is the fixed coordinate of the comoving test particle, and r is the 
proper time of the central observer. However, in some cases, the dependence of 
the relative velocities is exclusively through the Hubble speed parameter v, 


= v { t , x) = d(r)x. 


(55) 


Remark 4. The expression (55) for v is the Hubble speed of a comoving test 
particle with curvature-normalized coordinates (r, x) ■ However, it is important 
to recognize that the relative velocities that we calculate in this section, as func¬ 
tions of V, are those of test particles located at different spacetime points, of the 
form {t,x)y where t < t. 
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If there is a one-to-one correspondence between Hubble velocities to geometric 
velocities (Fermi, kinematic, spectroscopic, or astrometric), then an addition 
formula for the geometrically defined relative velocities may be derived from 


Eq.(54) as the following theorem shows. 


Theorem 5. Let Vg denote one of the four geometrically defined relative veloc¬ 
ities considered in this paper, and let v denote the Hubble velocity of a test par¬ 
ticle. Suppose further that Vs has a 1-1 functional dependence on v, Vg = fiv). 
Then the geometric velocity, Vs 3 , of (32 relative to (do o-t proper time Tq of /3o 
can be written in terms of Vgi and V 2 s, the geometric relative velocities of j3i 
relative to (3 q and P 2 relative to /3i, respectively. Here Vgi is the relative velocity 
of jdi at proper time Tq (it is possible for Tq = Tq ) of (3 q and Vg 2 is the relative 
velocity of /?2 at proper time ti of jdi. The relationship is 


, I a(T+) 1 a To+ 1 

Vs3 = f{TT^f {Vsl)+^V^f {Vs2) 
.a(To ) a[Ti) 


(56) 


Proof. The Hubble velocity, V 3 , of (32 relative to (3 at Tq can be expressed as. 


V3 = a{To)X2 = a(To+)xi -f d(Tp+)(x2 - Xi) 


(57) 


d(To+) a(ro+) 


= )Xi + 

a[To ) a[Ti) 

= 

a(To ) a{Ti) 


a(Ti)(X2 - Xi) 


where Vi is the Hubble velocity of /3i relative to Pq at proper time Tq of Pq and V 2 
is the Hubble velocity of P 2 relative to Pi at proper time ri of Pi. By assumption 
Vsi = f{vi) and Vs 2 = f(v 2 ). Since / is 1-1, we may write vi = /“^('Csi) and 
V 2 = f~^ivs 2 ). Substituting these expressions into Eq. (57) we obtain. 


j-i 

a{To )' 


1^3 = ^^r\vs 




d(ri) 


/ (^ 82 )- 


Then from Eq. (58), we have that. 


Vs3 = fiv3) = f 


a[To ) 


a(T-o^) .-1 


d(Ti) 


/ (^^^ 2 ) 


(58) 

(59) 

□ 


While Eq. (56) is less general than the other formulas developed so far, it is of¬ 


ten easier to use when working with a particular spacetime where it is applicable. 
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We conclude this section with a corollary that gives a kinematic velocity addition 
formula for Robertson-Walker spacetimes with scale factors following power 
laws, i.e., scale factors of the form, 


a{t)=t°‘, (60) 

where a > 0, a ^ It was shown in [3] and [5] that the proper radius 
of Mt is an increasing function of r and. 


_ PMAa) 

- r(^) 


(61) 


Let Fa{z) = 2 F 1 (|, ^;z^) where 0 < 0 < 1, and where 2 -f^i(-, S S 0 

is the Gauss hypergeometric function. Define Ga{v) = {F~^ {Ca + 

Specializing to the case x > 0 convenience only), for a comoving test particle 
(see 0) then, 

Wkin = fiv) = \/l - G2“(u), (62) 

where the function / plays the same role here as in Theorem For 0 < a < 1, 
V is bounded above by j^Ca, but for a > 1, x and v have no upper bounds 
[S]. It is readily seen that. 


V = f ^(Ukin) 



(63) 


Corollary 2. Let the scale factor for Robertson-Walker spacetime be given by 
a{t) = with a > 0, a 7 ^ 1. For the scenario depicted in Figure with 
X 2 > Xi > 0; ^^6 kinematic velocity addition formula is given by, 


VkinS = + (1 - 


(64) 


wherein we refer to Eg. (631. 


Proof. From (13) we obtain. 


a(ri) 


and thus. 


«(r„) = 


„a-l 


d(ro) ^ 


(65) 


( 66 ) 


The specialization of Eq. (561 to the kinematic relative velocity in the configu¬ 
ration depicted by Figure 1 is. 


®The case a = I gives the Milne universe 
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( 67 ) 

□ 


?^kin3 = / ( / ^(^^kinl) + ^(Wkin2) 

a(Ti) 


Combining this with Eqs.(62) and (66) then yields Eq. (64). 


We mention that the cases of negative % values, and/or xi > X 2 are also easily 
handled at the expense of some absolute values, and that similar results for the 
cases of Figures and may be obtained using the same method. 

Remark 5. The velocity addition formula given by Corollary is more spe¬ 
cialized than the formula given in Theorem [7| but, as in special relativity, has 
the feature that there is no explicit time coordinate dependence, unlike the more 
general velocity addition formula given in Theorem^ 


9 Examples 

In this section we apply the results of the previous sections to find explicit 
expressions for the velocity addition formulas of comoving observers in particular 
Robertson-Walker spacetimes. 

9.1 The de Sitter universe 

The Hubble velocity for a comoving test particle in the de Sitter universe is 
given by, 


V = v{t, x) = Hoe^‘>'^x- (68) 

where Hq is the Hubble constant. Expressed in the notation of Theorem 
it was shown in |1] that the kinematic velocity of a test particle relative to a 
central observer is given by 


Wkin = f{v) 


V 

Vi + v^’ 


(69) 


for |w(t, dl < V— 1 and |r'(-,x)| > =. Then, 

' ^ I V -A/I 


V = f (Ukin) = 


■Ckin 


\/l - 


(70) 


We use Theorem to find a velocity addition formula for the kinematic relative 
velocity in the scenarios depicted in Figs. [^and[^ Applying Eq. (57) to the 
scenario depicted in Fig. [l]we have that 


V3 = iJoe^“^“x2 = vi-\- 

a(ri) 


(71) 
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Using (131 we have, 


a{To) 


c )-^^0 ”^0 


0 (^ 0 ) 


a(Ti) e-f^oTi a{Ti) ’ 


( 72 ) 


Combining Eqs. and ( [T^ and substituting them into the formula Vkins = 

ds the addition formuh 

^kinl \/l ^kin2 ^kin2 


along with Eq. ( |70[ ) yields the addition formula. 


?^kin3 = 


1 + 2r)kinll^kin2\/l — ' 


^kin 2 


Turning to the scenario depicted in Eig. we have by Eq. (571 that 

a(t+)^ 


a(t+ 

^3 = ■ 


a{t ) a(ri) 


Vl. 


(73) 


(74) 


Again, using (13) we have, 

g(T- 2 ) 

a{t+) 

a{T2) 


aHoT2 




pHoTi 


Q('ri) ^ _ 

a{t~) e^o‘” 


II 

\/i-^L3 

(75) 

II 

V1 - ^kin2 

(76) 

a(Ti) _ 

a{t-) 


(77) 


Substituting Eqs. (751, ( [7^ , ( [77| into Eq. (74) and then using Eq. (70) gives 
the velocity addition formula. 


WkinS = "f^kinl ^ ~ ■^kin2 + ^kin2 • 


(78) 


9.2 Addition Formulas for a{t) = tl/3 

For this scale factor. 


X 

3t2/3 ’ 


and it was shown in [4] that. 


and 


Vkin = V 


'CFernii — 1^(1 “b U ), 


(79) 

(80) 

(81) 
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with |z;| < 1. Since the function defined by Eq. (811 is 1-1, the Hubble velocity 
can be expressed in terms of the Fermi relative velocity as 


'^(^Fermi) — \ «l^(^Fermi)5 

^V^Fermij ^ 


where 


C'(wFermi) = 


— 27?;Fermi + -\/27(4 -|- 27Upgj.jj,j) 


1/3 


(82) 


(83) 


We apply Theorem to the scenario depicted in Fig. to obtain addition 
formulas for the Fermi and kinematic relative velocities. From Eqs. (131, 


(S) 


1/3 


a(ri) 


/I-»Lr 


Combining this with Eq. (80) gives, 

n = (1 - 


Eq. (57) for this situation therefore reads, 

a(To) 


V3=Vi + -r^-^V2 = Hi -k (1 - Vl)v2. 
aiji) 


(84) 


(85) 


( 86 ) 


For the Fermi relative velocity it follows from Eq. (86) that 


^'Fermi3 = [^'(^'Fl) + (1 “ ^^(i’Fi)^)w(hf 2)] [l + (h(hfi) + (1 “ h(hfi )^)h(hf2 )) ^] , 

(87) 

where v{v-pj) = ^(HFermy) IS given by Eq. ([8^ for j = 1, 2. 

For the kinematic relative velocity, we obtain the velocity addition formula by 
use of Eq. (80|) in Eq. (1^. The result is. 


Hkin3 — Hkinl T (1 )Hkin2 ■ 


( 88 ) 


9.3 Addition Formulas for a{t) = ^1/2 

For the scale factor a{t) = for the radiation dominated universe, the Hubble 
speed is given by. 


V = 


X 

2 ^' 


(89) 
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and it was shown in [1], that the kinematic and Fermi relative velocities of a 
comoving test particle are given by, 


Wkin = sm V 


and 


i^Fermi = (cos z; + f sin v) sin z;, 

where |z;| < 7r/2. From Eqs. for the scenario depicted in Fig. 


fe) 


1/2 


_ /l 


a(To) ^kinl- 


Combining this with Eq. (90) gives, 


Therefore, 


Ti = To cos ZZl. 


, «(^o) 

+ T7- vV2 = Vi +V2 COSTi 

a(ri 


(90) 

(91) 

(92) 

(93) 

(94) 


For the kinematic relative velocity addition formula, we then have that 


■^kins = sinzza = sin(z;i + V 2 coszzi) 

= sin(z;i) cos(z ;2 coszzi) + cos(z;i) sin(z ;2 coszzi) 

= zzkini cos (^sin“^(z;kin 2 )(95) 
+ \/l-^LiSin (sin-^(z;kin2)Y^l-^^Li) • 


For the Fermi relative velocity, since the function defined by Eq. (91) is 1-1 on 
the interval (—7r/2,7r/2), in principle an inverse function can be found that gives 
the Hubble velocity as a function of the Fermi relative velocity. By Theorem 
we can then find an addition formula for the scale factor a{t) = 


Using the results of 015], general velocity addition formulas can be similarly 
formulated for scale factors of the form, a{t) = for a > 0. 


10 Conclusion 

Natural generalizations of the the special relativistic velocity addition formula 
arise from the consideration of comoving observers and particles in Robertson- 
Walker spacetimes. As described in the introduction, such a configuration is 
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the natural generalization of special relativistic inertial frames in standard con¬ 
figuration, related by Lorentz boosts. 

The introduction of geometrically defined relative velocities, Fermi, kinematic, 
spectroscopic, and astrometric came about following discussions about the need 
for a strict definition of “radial velocity” at the General Assembly of the Inter¬ 
national Astronomical Union (lAU), held in 2000 (see, e.g., Given 

the observation of the preceding paragraph, it is natural to investigate addition 
formulas for these geometric relative velocities. 

General velocity addition formulas for the Fermi and kinematic relative veloc¬ 
ities, associated with spacelike simultaneity were found in Section and the 
analogs for the spectroscopic and astrometric relative velocities associated with 
lightlike simultaneity were developed in Section In this degree of generality, 
with the important exception of the spectroscopic relative velocity (see The¬ 
orems and 1^ , the velocity addition formulas depend, not only on the space 
coordinate x (in curvature coordinates) of the comoving observers and particles, 
but also on their time coordinates. 

However, if a geometrically defined relative velocities is a bijective function of 
the Hubble velocity, in the sense described in Section]^ then the dependence 
on time coordinates in the addition formulas can, at least in some cases, be 
eliminated. This is carried out for the de Sitter universe and general power 
law cosmologies with scale factors of the form a{t) = t“,a > 0 (a = I gives 
the Milne universe, i.e., special relativity for which the usual velocity addition 
formula is recovered by our methods). Simple examples, including the radiation 
dominated universe, are given in Section 

When a > 1, the Robertson-Walker cosmologies with scale factor a{t) = in¬ 
clude event horizons. Fermi and optical coordinates cannot be extended beyond 
the event horizon [5] , so the geometric relative velocities of test particles in that 
region of spacetime are undefined. By contrast, Hubble velocity fields beyond 
the event horizon have been studied m and are of interest as part of the large 
scale structure of the universe. A possible future direction for research could 
be to investigate whether the velocity addition laws in this paper could be used 
to extend the definitions of geometric relative velocities through the use of an 
intermediary observer, within the central observers event horizon, in a way that 
is well defined, and for more realistic spacetimes. 

Acknowledgment. J. Reschke was partially supported during the course of 
this research by the Interdisciplinary Research Institute for the Sciences at Cal¬ 
ifornia State University, Northridge. 
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